Degeneracies in the length spectra of metric graphs 

U. Gavish^ and U. Smilansky^ 

^1 1. Laboratoire de Physique Theorique et Hautes Energies, CNRS UMR 7589, 

("^ ■ Universite Paris 6, 4 Place Jussieu, 75005 Paris, France 

CN I 2. Weizmann Institute of Science, Rehovot 76100, Israel 

Abstract. The spectral theory of quantum graphs is related via an exact trace 
QQ . formula with the spectrum of the lengths of periodic orbits (cycles) on the graphs. 

The latter is a degenerate spectrum, and understanding its structure (i.e., finding 

out how many different lengths exist for periodic orbits with a given period and the 

rS^ , average number of periodic orbits with the same length) is necessary for the systematic 

'T^' study of spectral fluctuations using the trace formula. This is a combinatorial problem 

<^ . which we solve exactly for complete (fully connected) graphs with arbitrary number of 

j^ I vertices. 
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1. Introduction 

The interest in the spectral properties of the Schrodinger operator on metric graphs 
(known also as "quantum graphs") increased dramatically after it was found that 
quantum graphs provide an excellent paradigm for the study of spectral fluctuations 
in quantum chaotic systems [21 |3]: The spectral density of quantum graphs can be 
expressed as an exact trace formula [H [2] in terms of the spectrum of the lengths of 
its periodic orbits (PO) (also called cycles) which is analogous to the asymptotic semi- 
classical trace formula |4] . Moreover, the quantum graphs posses a Liouvillian analogue, 
which, under some well understood conditions is ergodic. At the same time, extensive 
numerical simulations and tests can be performed with rather modest computational 
effort allowing detailed comparison of the spectral statistics with the prediction of 
Random Matrix Theory, and study the systematic deviation from it. The simple finite 
graphs are essentially one dimensional (albeit not simply connected) systems. They 
display spectral complexity under one important condition: the lengths of the bonds 
must be rationally independent. 

The main tool in the theoretical discussion of the spectral statistics of quantum 
graphs is the above-mentioned trace formula. It can be written explicitly as 

d{k) = E ^(^ - ^.) = 7 + E E 4"^e^''^ ' (1) 

where Ej = k'j are the eigenvalues of the Schrodinger operator. The kj's form the 
wave-number spectrum. C is the total bond length, i.e., for a graph with B bonds, of 
lengths Lb, it is given by C = Y.b=i Lb- Vn denotes the set of PO's of period n. Each 
periodic orbit contributes a term which consists of a "transition amplitude" A^J^\ and a 
unimodular factor with a phase which is determined by the length of the corresponding 
n-bond PO 

B B 

lp = ^qbLb ; gt e {0,1,2, ■■■} and E * = '^ • (2) 

b=l b=l 

The length spectrum is highly degenerate. Each degeneracy class contains all orbits 
which traverse the same bonds the same number of times, but not in the same order 
(up to cyclic permutations) and only these (since the bond lengths Lb are rationally 
independent). That is, a degeneracy class of n-bond PO's consists of orbits which have 
the same code q*^""^ = (gi, ■ ■ ■ , (js), with Y,b % = n. Not every set of nonnegative integers 
{qb\ the sum of which is n represents a degeneracy class - the graph connectivity and 
the periodicity restrict the possible codes. The trace formula can be written as 



^ n „(n) 



-^'iklf^ 



(3) 



_peq(") 

Where the contributions of the different orbits in the degeneracy class q*^"-* were lumped 
together to the sum in the square brackets, all having the same phase factor. 

There are two prominent examples where a detailed information about the 
degeneracy classes is needed. The first example emerges in attempts to understand 
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the conditions under which the spectral fluctuation of a quantum graph follow the 
predictions of Random Matrix Theory. A standard tool is the computation of the 
spectral autocorrelation function 

where d{k) = d{k) — - is the fluctuating part of the spectral density and the domain of 
integration [k ~ A,k + A] is arbitrarily large. Substituting the explicit expression ([3]) in 
(jl]) one sees that the autocorrelation function depends on the squares of the individual 
contributions of the degeneracy classes (the terms in square brackets in ([3])). 

The second example is encountered in the context of " hearing the shape of a graph" , 
that is, in attempts to reconstruct the connectivity and the length spectrum from the 
energy eigenvalue spectrum of the quantum graph [9]. The main tool is again the 
trace formula ([3]) or rather its Fourier transform d{l) = J dkd{k) exp{ikl) . d{l) is a 
distribution supported on the length spectrum, with weights which can be read off from 
([3]). The length spectrum (its composition and weights) is therefore useful to obtain the 
information about d{l) and in turn also about the connectivity of the graph. 

The leading asymptotic (for large n) contribution to the number of degeneracy 
classes in a general connected graph was obtained by Berkolaiko [5J. The number of 
degeneracy classes and the number of PO's in each class were obtained by Tanner [6] 
for binary graphs up to order 6. 

In this paper we present an exact expression for the number of classes for fully 
connected (complete) graphs of any order. 

We start by defining precisely graphs, PO's, and their degeneracy classes. We then 
compute the number of degeneracy classes for general fully connected graphs, the total 
number of n-bond PO's, and obtain the mean degeneracy of the classes as the ratio 
between the two. Finally, we present numerical results and interpret them. 

2. Graphs, Periodic orbits, Degeneracy classes 

A graph, G, of order V, is a set of V numbered vertices, some of which are connected 
by a bond (no more than one bond between two vertices, no bond connects a vertex 
to itself). The number of bonds connected to a vertex is the vertex valency. The 
connectivity matrix of G is defined by 

„ ,^, r 1 if i ciiT'd j are connected ,^. 

■'^ ^ ^0 otherwise. ^ ^ 

A fully connected simple graph Ky is a graph where each vertex is connected by a 
single bond to any other vertex (beside itself): Gij = 1 — 6ij, i,j G G. 

A periodic orbit (PO) on a graph G is a sequence of vertices, [fi,f2, ...,fn] with 
Cvi,v,+i = 1 and Vi = Vn- PO's that can be obtained from one another by a cyclic 
permutation of their vertices will be considered identical. 

Consider an integer set {qh}h=(ij)={j,i), i,j G G with J^Qb = n. A degeneracy 
class of n-bond periodic orbits is a set of all the n-bond PO's each of which passes 
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1^4^2^1^4^5^1 

2 1^4^2^1^5^4^1 

1^5^4^2^1^4^1 

1^2^4^5^1^4^1 



1.d 








1^2^4^1^4^5^1 



y' y y y y -^ y 



Figure 1. The fully connected graph K^ and one of its degeneracy class 

exactly qb times over the bond b. All these PO's are of the same length and since 
the bond lengths are rationally independent - all PO's of the same length belong 
to one class. The degeneracy of a class is the number of distinct PO's in it. Fig. 
l.a shows the fully connected graph G = K^. Fig. l.b shows the degeneracy class 
(q'(i,2) = (1(2,4) = ?(4,5) = 1(5,1) = 15^(1,4) — 2)- Fig. l.c lists all the PO's in this class and 
Fig l.d shows two of them explicitly. 

Let Nc{n, G) be the number of classes of n-bond PO's in G and Np{n, G) the total 
number of n-bond PO's in G. The mean degeneracy of n-bond PO in G, Dn{G), is 
defined by: 

In the next section we shall provide exact expressions for Nc{n,Kv), Np{n,Kv), 
and Dn{Kv), that is, for the number of classes and PO's, and the mean degeneracy of 
fully connected simple graphs. 

3. The number of classes of n-bond periodic orbits on a fully connected 
graph with V vertices 

Since a fully-connected graph is determined uniquely by specifying V, we may use more 
compact notations, replacing Dn{Kv) -^ D{n,V) and so on. We start by obtaining 
Nc{n,V). 

Let V <V and Nc^^{n,v) be the number of n-bond degeneracy classes in K^ which 
contain PO's that passes through all the v vertices. Note that Nc,v{n,v) < Nc{n,v) 



2.a 







Figure 2. 4-bond degeneracy classes of the fully connected graph K4 



because not all PO's pass through all the vertices. More precisely, 



Nc{n, V) = Y: 



\Nc,v{n,v) 



(7) 



v=l 



that is, the number of classes is a sum of the number of classes with PO's that use 
exactly v vertices. The factor (j accounts for the possibilities of choosing these v 
vertices. All such choices have identical contribution since Ky is fully connected. 

For example consider n = V = 4. Figs. 2. a, 2.b and 2.c show all the 4-bond classes 
of Ki grouped according to the sum in Eq.(I7j). There are (2) = 6 ways for choosing 
2 vertices (see Fig.2.a), f j = 4 for choosing 3 vertices (corresponding to each line 

in Fig.2.b), Lj = 1 ways for choosing 4 vertices (Fig.2.c). These figures shows that, 
iV,,i(4,l) = 0, iV,,2(4,2) = 1, iV,,3(4,3) = 3, iV,,4(4,4) = 4. 

The reason why we have expressed Nc{n,V) in Eq.Q in terms of Nc^^{n,v) is 
because the latter was calculated by R. C. Read in Ref. ^0\. For completeness, a 
concise derivation along the lines of that work is given in Appendix 1. The result 
derived in Appendix 1 is that Nc^^{n,v) is v\ times the coefficient of x^t"' in the Taylor 
expansion (near t = x = 0) of In {E{x, t)), 

■ln(E(x,t))U=o,t=o, (8) 



Nr„ 



n,f 



n\ dx" df^ 



where 



X 



E(x,t) = E2-'^-(l-t) 



^-\v{v-l) 



v=Q 




(9) 
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From this and Eq.([7]) one sees that the number of classes of n-bond periodic orbits on 
a fully connected graph with V vertices is 

An explicit expansion of E{x,t) yields 

E{x,t)= Y. En^vt-'X^ (11) 



n,i'=0..oo 



with 



(-1)^ (v\ (f, + s{v-s)~iyn~f^+ (I) + (--) - 1\ 

(12) 



,f^,^ 2^t;! \sj\ /i J\ n- ^ 



fi=0..n 

One can also derive the following recursion relation (see Appendix 2) 

(v - 1)! 

Nc,v{n,v) =v\En,v- Yl Yl /, _ '^ Nc,v{m, k)En-m,v-k (13) 

which enables a fast numerical calculation of Nc{n, V). 

As an example, consider the case V = n = 4. Substituting these values in Eq. lTTO!) 
one obtains Nc{4:, 4) = 21. This result is confirmed in Figs.2.a-c which present all the 21 
classes of PO's. Eq. A.l, in Appendix 1 of Ref. [5] provides the asymptotic behavior, 
for large n, of the number of classes N^. (for any connected graph): 

nB-V+l B~l 1 

N,{2n, V) + N,{2n + 1, y) ~ -^^—^(1 + 0(-)), (14) 

where B is the number of bonds, for the complete-graph one has B = V{V — l)/2. By 
substituting this value into Eq. flT^ we verified numerically that the asymptotic behavior 
of Eq. (fTO|) for n ^ V matches Eq. (fT^ . The results are shown in Figure 3 which presents 
the ratio {N^{n, V) + N^{n + 1, V))/{N^'y"'P{n, V) + N^'y'^P{n + 1, 1/)), for even n. The 
superscript asympt stands for the values obtained using Eq.([T] 



4. The mean degeneracy of n-bond periodic orbits on a fully connected 
graph with V vertices. Numerical results. 

To obtain the mean degeneracy Eq. ([6]), we need also to derive an expression for Np{n, V) 
i.e. the number of n-bond PO's in Ky- Let us first calculate the number of n-bond closed 
trajectories. The number of n-bond closed trajectories is given by 

N{n, V) = TrC" (15) 

where C is the connectivity matrix defined in Eq.(l5]). In our case it is given by: 

a,,{Kv) = 1 - <5.„ (16) 

and its eigenvalues are: Ai = V^ — 1, A2 = A3 = ...Xy = —1- From Eqs. flTSl) and flT6|) it 
follows that 

N{n,V) = {V-ir + {V-l){-ir. (17) 
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N^{n,V) + N^{n + XV) 
N';'"""{n,V) + N';""'"(n+^,V) 




Figure 3. Asymptotic behavior of the number of degeneracy classes 

N{n, V) is the number of ra-bond PO's but is different from Np{n, V) in Eq. ([6]) since 
in the latter, PO's that can be obtained from one another through a cychc permutation 
are considered to be the same PO. 

For simphcity, let us assume that n is a prime number, thus avoiding the 
complications arising from the presence of PO's which are repetitions of a shorter PO. 
With this assumption, each of the PO's counted in Eq. flT7|) is one of n PO's that can 
be obtained from one another by cyclic permutations of the vertices. Since we regard 
all such cyclicly-equivalent PO's to be the same one, Np{n, V) and N{n, V) are related 
by 

Np{n,V) = -N{n,V). (18) 

n 

From Eqs. ([6]) and (TTSi) one has 

c(„,y) = '";i)y^.-f + '-1)""^-^)) (19) 

J2v=l [v)d^dt^^^(^(^^'^)) \x=0,t=0 

The mean-degeneracy D{n, V) and its logarithm are shown in Figures 4, 5 and 6. These 
plots where generated using either Eq. (fT9l) or the recursive relation Eq. (fT3l) . Figure 
4 presents the n-dependence of D{n,V) for fixed values of V. Figure 5 shows the In- 
dependence of D{n,V) for fixed values of n. As seen from these figures, the mean 
degeneracy grows rapidly and achieves values much larger than 2 already for small 
(i.e., much smaller than the number of bonds) values of n. On the other hand, in the 
limit V/n ^ oo it approaches 2 (most classes contain only a single PO and its time 
reversal). Approximating D{n, V) ~ N}2nvt+N(2n+iv) ^^*^ then using the asymptotic 
expression Eq. flT^ together with Eq. flTTl) one has 



D{n, V) ^ V{V' -V- m^'-'^-j^ (20) 



n 



y(v-i)/2 
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Mean Degeneracy at Fixed V 




9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 



Figure 4. The mean degeneracy for fixed number of vertices, V. 



Taking the logarithm of both sides and keeping only terms containing n (assuming 
logn -C V) one obtains 



\og{D{n, V)) ^n\og(y - 1) - -ViV - 1) logn. 



(21) 



taking the derivative with respect to V yields the approximate value of V in which the 
maximal mean-degeneracy is obtained: 



\ rn 



n 



logn 



(22) 



Although this estimation was derived for large n it shows reasonable agreement with the 
peaks in Figure 6. For example, the two maxima marked with arrows, for n = 20 and 
n = 30, are located in the vicinity oi V = 3.9 and V = 4.5 respectively, in agreement 
with Eq.( 
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Appendix 1. Derivation of Eq. ^ 

In this appendix we derive Eq. ([H]) along the lines of Ref. [TD]. 
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Mean Degeneracy at Fixed n 




15 ?? 29 36 43 50 57 64 71 
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Figure 5. The mean degeneracy at fixed PO's lengths. 
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Figure 6. The logarithm of the mean degeneracy at fixed PO's lengths. 



Definitions A graph of order V is a set of V numbered vertices some of which are 
connected by a bond (no more than one bond between two vertices, no bond connects a 
vertex to itself i.e. no loops). The number of bonds connected to a vertex is the vertex 
valency. Fig. 7.1 shows a graph of order 6. 

A multigraph is similar to a graph except that there can be more than one bond 
between two vertices (Fig. 7.2). A graph is a specific case of a multigraph. Let q{ij){g) 
be the number of bonds connecting the vertices i and j in a multigraph g. We shall refer 
to q(ij){g) as the degree of the bond {i,j) in g. For example, in Fig. 7.2 ^(2,3) = ^(3,2) = 3. 
Note that the pairs (i,j) are not- directed, that is, (i,j) and (j,i) are considered as the same 
bond. 
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7.2 A multigraph 




7.3 A connected multigraph 





yy A Multigraph not 
^ embedded in Fig. 
7.3. {Here: q5 2='' - 
in Fig. 7.3 ^^^=0) 




7.4 Fully-connected graph 



7.5 A Multigraph embedded in 7.3 



7.6 



7.7 



1 ^2 ^4 ^1 

":- ■ ■ V- ■ \>- 

2^4 ^1 ^2 





7.8 An even multigraph 



7.9 A connected even multigraph 



Figure 7. Graphs and Multigraphs 



A multigraph is connected if, by moving on the bonds, one can pass between any 
two of its vertices (in particular, all valencies are > 1, Fig. 7.3). A graph, Ky, is fully- 
connected if each of its V vertices is connected to all other vertices (Fig. 7.4). Thus, Ky 
has V{V - l)/2 bonds. 

Let g he a multigraph of order v. g is said to be embedded in a graph G if it can be 
obtained from G by first adding and deleting bonds between vertices that are connected 
in G, and then deleting some of the vertices which have zero valency (now, after the 
addition and deletion of the bonds). By this definition, the order of G is larger or equal 
to V. If q{i,j){G) = 1 then q{ij){g) = 0, 1,2... and if q{i,j){G) = then q(i,j){g) = 0,. The 
multigraph in Fig. 7.5 is embedded in the graph of Fig. 7.3 while the multigraph in 
Fig. 7.6 is not. 

A trajectory is a sequence of vertices, adjacent pairs of which are connected. If it 
is closed, i.e. it starts and ends at the same vertex, the trajectory is a PO. Actually, 
one can associate several PO's with each closed trajectory since one can start in any of 
the trajectory points, however, we shall refer to all of these as a single PO that is say 
for example that 1^2^4^1is the same PO as2^4^1^2 and two PO's are 
distinct only if they can not be obtained from one another by such a cyclic permutation. 
Thus, the two closed trajectories in Fig. 7.7 are the same PO. 

An even multigraph is a multigraph where each vertex has an even valency (Fig. 7. 8). 
For any connected even multigraph (often called Euler multigraph) one can always find 
a PO that passes on each bond exactly once (Eulerian circuit). Often, there is more 
than one. For example, the two PO's (Fig.l.d) 1^2^4^1^5^4^1 and 
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1^2-^4— >l—i>4—i>5—i>l passes once on each bond in the connected even 
multigraph Cg shown in Fig. l.b. 

A class of n-bond periodic orbits, Cn, in a graph G is an n-bond connected even 
multigraph which is embedded in a labeled graph G. Gn, is specified by specifying the 
set of bond-degrees {q(^ij){Gn)}i,jeG, where J2i,j j>iQ{i,j){Cn) = n. A PO is said to be 
in Gn if it consists of n steps passing exactly q(^ij){Gn) times between i and j. By this 
definition, all PO's in C„ have exactly the same length, independently of the choice of 
bond lengths and therefore, at a given energy, have the same action. Fig. l.b shows a 
class of 6-bond PO's, Gq, embedded in the fully-connected graph G in Fig.l.a. The PO's 
in this class are listed in Fig. I.e. Two of them are drawn in Figs.l.d. 



Proof of Eq. ^ The proof of Eq. ([8]) is based on that given in Ref. [TO] which 
can also be used to treat the case of classes of graphs with and without loops, and 
multigraphs with loops. Consider a set of n labeled vertices, and the set f2(n, v) = {Gn,v} 
of multigraphs with n bonds that one can draw on these v vertices {on means using all 
of them - i.e. these multigraphs are of order v). To each of the v vertices we assign a 
sign, +1 or -1. There are 2^ such possible assignments. For a given assignment S, we 
define the sign of each bond in Gn,v to be the product of signs of its two vertices. The 
sign, a{Gn,v, S), of a multigraph Gn,v £ fl{n,v) is then defined as the product of signs 
of all its bonds. Thus, 

CT{Gn,v, S) = (_l)^-(Gn,.,5) _ (_i)MGn,.,5) (23) 

where V- is the sum of valencies of the negative vertices and /i the number of negative 
bonds. The sum of the signs of Gn,v for all possible S is Y.si~'^)^~^^""''^- Summing this 
over all members oi Q{n,v) one has 

Y^ [ y^(_X)V'-(G„,.,5)^ y- / y- ^_iY((Gn,v,s))] (24) 

G„,ve^(n,v) \ S J S \G„,ven{n,v) J 

In the right hand side the order of summation was reversed and Eq. (1231) was used. 
Consider the left hand side of Eq. fl2^ . If Gn,v is an even multigraph, then V_{Gn,v, S) 
is an even number for any S and therefore ^^(— 1)^-('^"''"'^) = 2'". If Gn,v is not even, 
then at least one of its vertices, say A, has an odd valency. Since for each assignment S 
in which A is negative there exists S' which is identical to S except that A is positive in 
it, and since a{Gn,v, S) = —a(Gn,v, S'), one has ^^(— 1)^-('^"."''^) = for any Gn,v which 
is not even. Thus, the left hand side of Eq. fl2^ is the number of even multigraphs in 
f2(n, f) times 2^. To obtain the right-hand side, consider the ( j assignments in which 
exactly s of the vertices are positive. The number of ways to put /i identical balls in 
s{v — s) identical boxes each of which may contain any number of balls, is (f^+^'^'"~^>~ j 
and therefore this is the number of ways the fi{Gn,v, S) bonds which join the s positive 
with the V — s negative vertices can be placed. The remaining n — fi bonds may be 
placed between the (2) + (^2'*) pairs of vertices with identical signs, that is in 

V n-fi J 
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different ways. (To enable compact writing, here and below, we assume that the binomial 
coefficients have the properties [Vj = for b > a and b ^ 0, and (Vj = 1 for any a). 
Summing over all possible ii{Gn,v, S) one gets the total contribution of all assignments 
in which exactly s vertices are positive: 

|:(-i)- (" + ■''" - ■') - 1) (" - " + 0_+ ("i") - M . (26) 

This contribution is the coefficient oft" in (^l—t)^^^/'^'>'"^'"~^'>(^j^)'^^'"~^\ Thus, the number 
of n-bond even multigraphs one can draw on v labeled vertices is given by v\ times the 
coefficient of t^-x"" in the power expansion of: 

We are interested in Nc^y{n,v), i.e. the number n-bond connected even multigraphs one 
can draw on v labeled vertices. It is a known result in graph enumeration theory that 
the generating function of the connected set of (labeled) graphs is given by the log of 
that of the non-connected set [11]. Thus, Nc^v{n,v), is v\ times the coefficient of t"'a:'' in 
the power expansion of lia{E{x,t)) which proves Eq.(IH]). 



Appendix 2. Proof of Eq. (E 
Define the expansions 

oo 

E{x,t) = Y.E,{t)x\ {2i 



v=0 

and 

CX) 



\n{E{x,t)) = Y.L,{t)^. (29) 



1 v\ 

v=l 



Nc,v{n,v) is v\ times the coefficient of x^t"- in ln{E{x,t)) and therefore 

oo 

L,it) = Y.N,,,{n,v)r. (30) 



n=0 



There exists a useful recursion relation between the coefficients in Eqs. ( [28]) and fl29l) : 
L,{t) = v\E,{t) - Y, |!^~||; Lfc(t)E,_fc(t) V > 1 

k=i yi^ ~ ^)- 

L,{t) = E,{t) (31) 

Expansion of Eqs. (12^ and (1501) in powers of t yields Eq. (IT^ . 
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